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Abstract 

The cross section for the reaction dd ^Herj cfose to threshold is 
estimated in a two-step model, similar to that which successfully describes 
near-threshold pd — > ^Herj production. A vr-meson, created in a dp ^ 
'^He 7r° process on one nucleon in the target deuteron, is converted into an 
ry-meson by a secondary 7r°n — > n?? reaction on the other nucleon. The 
kinematics are such that both processes are nearly physical so that only 
small Fermi momenta are required. The predictions of the model are in 
good agreement with recent Saclay experimental data. 

An analogous model, when applied to the near-threshold n ^He — > ^He rj 
reaction, predicts a cross section of a rather similar size, but experimental 
difficulties render this a doubtful means of studying the rj ^He interaction. 
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1 Introduction 



The strong energy variation observed in the pd — >■ ^He t] amphtude near threshold 
[|I], 3 lead to the speculation of the existence of quasi-bound ?7-nucleus states 
for much lighter nuclei than originally suggested [^]. This incited groups working 
at the SPESIV p and SPESIII ^ spectrometers at Saturne to investigate the 
next heavier nucleus through measurements of the dd ^Herj reaction. First 
results indicate a '^Her] threshold production cross section which is about a factor 
of 40 less than for ^Her], combined with a weaker energy dependence. 

The near-threshold energy variation is dominated by final-state-interaction 
effects, which reflect the ?7-nucleus dynamics, and these should be largely in- 
dependent of the particular entrance channel. In contrast, the prime aim of the 
present paper is to try to understand the factor of 40 between the two production 
rates, and this leads us rather to studying the reaction mechanism itself. 

Near-threshold production of a meson as massive as the t] necessarily involves 
large momentum transfers and models which rely on spectator nucleons yield 
much too low cross sections 0, H. In a previous paper hereinafter referred to 
as I, we presented and evaluated a two-step model for the pd — >• ^He rj reaction. 



based on ideas put forward by Germond [10 and Kilian and Nann The 



basic assumption is that a virtual pion beam is created by the interaction of 
beam protons with one of the nucleons in the target and it is this pion beam 
which produces the rj, through an interaction with the second nucleon in the 
target. A similar model has been independently studied in reference [0. 

The importance of this two-step mechanism is enhanced by what has been 
called the kinematic miracle by Kilian and Nann[|ll|]. They noticed that the 
kinematics are such that both sub-processes can proceed at very small Fermi 
momenta, with the constituent particles being almost on their energy shells. Any 
quantitative evaluation of such a model should not therefore be plagued by off- 
shell uncertainties and reliance on tails of Fermi momentum distributions. Taking 
advantage of the nearness to the physical domain, we expanded the intermediate 
state propagator to first order in the Fermi momenta to obtain a transparent 
physical estimate of the pd — > ^He t] cross section in terms of those for pp dir^ 
and 7r~^n — *■ r/p times simply evaluated form factors. 

Close to threshold the only significant energy dependence comes from the final 
state interactions. After taking these crudely into account, the model predicts 
cross sections which are too low by about a factoiQ of 2.8, which we attribute 
to the neglect of contributions of intermediate NN singlet states (d*). This is 
consistent with the observation that, with the same normalisation factor, the 
gross features of the threshold production rates for a;, rj' and mesons in the 
pd ^HeX reaction are equally well explained by the model [0. 

In the present paper, we extend the model by considering the dd ^Herj 

•^Due to a numerical approximation, the original factor of 2.5 in I should be corrected to 2.8. 
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reaction to be the result of a dp — > ^He process on one target nucleon, followed 
by 7r°n nrj reaction on the other. The predictions of this model in the forward 
direction near threshold, including final-statc-intcraction effects, are even closer 
to experimental data than for the lighter isotope, giving extra credence to the 
virtual pion beam concept. 

The formalism is completely analogous to that in I, so that the development 
given here can be brief. After presenting the generic form of the two-step diagram 
in §2, the amplitude is written as an integral over the ^He and deuteron Fermi 
momenta. By linearising the propagator in §3, the six-dimensional integral is 
reduced to a single one and an explicit prediction obtained for the dd — > '^Herj 
cross section. The effective wave technique used in I for handling the r^-nucleus 
final state interactions is generalised in §4 to the case of '^Herj. 

An alternative way of obtaining the same exit channel, in order to study the 
'^Her] final state interaction, is through the measurement of the n^He — > '^Her] 
reaction. This can be estimated in a similar two-step mechanism, and the method 
for doing this is sketched in §5, where an explicit formula is presented for this 
cross section in terms of the same pd — > ^He tt" and vr^n nr) cross sections that 
arise in the dd case. 

Our results and conclusions are presented in §6, where it is seen that the 
overall magnitude of our dd — > "^Her] predictions agrees even better with the 
experimental data than we found for pd — > ^He rj, possibly because of the absence 
of any Tinquantifiable d* contribTitions in the former case. On the other hand our 
treatment of the final state interactions leads to too sharp a threshold behaviour. 

Our n ^He —>■ ^He r] cross section estimate is of the same order of magnitude 
as for the deuteron-initiated reaction. Taking also into account the difficulty of 
working with a neutron beam or a tritium target, it suggests that this reaction 
is not competitive with dd — > '^He rj for producing low energy 77-mesons. 
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2 Nuclear matrix element 



In terms of nuclear matrix elements Ai with the same normalisation as in I, the 
cm. differential cross section for the unpolarised dd ^Her] reaction is written 

_ 1 P^^Y.\M\\ (1) 



dVt 16(27r)2sprf9 



spins 



with y/s being the total cm. energy of the rj ^ife-system and and pd the mo- 
menta of the 7] and deuteron in this system. 

There are strong restrictions on the amplitude structure at threshold for the 
rja system has there = 0~, which necessitates the deuterons to be in a L = 1, 
S = 1 state. The unique amplitude may be written in invariant form as 

1 



M = j^;^ V^m~m^ [iF (q' x ■ pa] , (2) 

where the q and q ' are the standard orthonormal sets of deuteron polarisation 
vectors. 

Though actually symmetric under the exchange of the two deuterons, at 
threshold we did not need to impose this in order to derive the spin structure of 
Ai in eq.(0). However in the forward direction, even away from threshold, it is in 
fact the only amplitude with the correct Bose symmetry. This means that in the 
forward direction, only deuterons with helicity m = ±1 can initiate the reaction 
and so there is a unique deuteron tensor analysing power, ^20 = +l/v^- This 
result is used in the analysis of the SPESIII experiment 0, which employed a 
tensor polarised deuteron beam. 

In terms of the amplitude F, the differential cross section of eq.(|l|) becomes 



18 (27r)8s Pd 



In analogy with the model proposed in I, one of the two-step contributions to 
the dd — > ^He t] is where a secondary vr" beam is created through the dp —>■ ^He 
reaction on the proton in the target deuteron. The tt" is subsequently transformed 
into an 7] by a ir^n —>■ rjirt interaction on the target neutron. This model is 
expressed in diagrammatic terms in fig. 1. 

The nuclear matrix element for the reaction can be written as 

M = NrV^^d f^d'k [^d^-^ ; ^. Mk, (4) 

J Eh J En [2txY 2E^ Eq - Eint + It 

where Eq = E^ + E^ is the total energy of the initial (and final) state. En and 
Eh are the energies of the neutron and helion in the intermediate state, where 
the total energy is Eint- 



4 



A^x is the product of individual matrix elements; 
Mk = Y.^ (^'^c ^ «/ )M (7r°n ^ n'r^j) M {di p^n) M {d'iP, n^K) , (5) 

int 

where the sum runs over the spin degrees of freedom of the intermediate particles. 
An index c signifies that the particle is represented by its charge conjugate spinor. 

In addition to the process described by eq.(|^), there is a contribution of similar 
size coming from the graph where protons and neutrons are interchanged, and 
two more contributions from charged intermediate pions. The latter are larger 
by a factor of two and so the amplitude of fig. 1 must be increased by an isospin 
factor of six. However, there are in addition an identical set of graphs where 
beam and target deuterons are interchanged and, in the forward direction, these 
ensure the correct m = ±1 helicity selection rule even away from threshold. Thus 
the m = ibl amplitudes of eq.(^) must be multiplied by a total combinatorial 
factor of Nt = 12, whereas the m = amplitudes vanish. 

We now study separately the four matrix elements of eq.(^) 



2.1 The deuteron wave function 

The momentum space wave function for an initial state deuteron at rest is 



M{di -^pcu) = T] 



V2 



V • (6) 



Here ?7„ is the neutron spinor and T]pc the charge-conjugate proton spinor, 
and ^sio) and ^d{(1) the deuteron S- and D-state components in momentum 
space. When the deuteron is moving with energy Ed-, the longitudinal part of the 
argument of the wave function is changed to 



(7) 



where 7 = Ed/rrid. 



2.2 The '^He wave function 

In the reaction model of fig. 1, we require the projection of the '^He wave function 
onto a bound state of a neutron and a ^He{h). In terms of the neutron and 
charge-conjugate helion spinors, this component of the final ^ife wave function in 
momentum space is 

M{nK - «/) = '^^/KClraik) Vn , (8) 
where k is the relative neutron-helion momentum inside the ^He. 
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If we take a completely symmetric Gaussian four-particle wave function for 
the '^He nucleus, then the relative wave function between the neutron and helion 
is also Gaussian, 

Mk) = ( ^ j exp(-|/?2p) . (9) 

The parameter j3 may be fixed by the point rms radius of '^He, 

< (1.46)2 = . (10) 

This gives f3 = 1.59 fm. 

In this simple model, the total number of p and n ^He pairs in the a-particle 
is No, = 4. 



2.3 Pion production pd 



7T 



3//e 



There are only two independent pd 7r° ^He amplitudes in the forward/backward 
direction |l^ . 



M{pcd ^ n'^K) = r]l^ [A {ea-pd) -iBa ■ (q x p^)] (he > 



(11) 



where pa is the direction of the deuteron cm. momentum. 

This decomposition is invariant in all systems connected by a boost along the 
beam direction. A represents the amplitude for helicity m = deuterons and B 
that for m = ±1 such that the aligned pd n^^He cm. cross sections are given 
by 



^ipd-.rr''He) 



mpfUh p^ |^|2 



4 (27r)2spd pp 

mpTUh IK 
4 (27r)2spd Pp 



B\ 



[m 



(m 



0) 



±1^ 



(12) 



2.4 (tt, 77) transmutation; ir^n rjn 
The unique threshold amplitude for n^n —>■ rjn is 



1 



VP J Vn 1 



(13) 



where s^p is the square of the r^-proton cm. energy. The amplitude for the more 
measurable vr^p rjn reaction is a factor of — larger, and the associated 
unpolarised production cross section is given by 



niz 



4 (27r)2s^„ p^ 



^ I ^ ( Srin ) 



(14) 
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3 Linearisation of the pr9pagator and 
evaluation of the matrix element 

Just as for the pd ^Herj reaction discussed in I, our two-step mechanism for 
dd ^Her] is such that the intermediate state is almost on the energy shell, 
corresponding to a nearly vanishing denominator in eq.@, even for small Fermi 
momenta q and k. A good approximation for the integral may then be obtained 
by expanding this denominator around the point g* = 0, A; = 0, retaining only 



linear terms. All dependence of the pd — > vr He and 



rjn amplitudes on 



the Fermi momenta are neglected when evaluating the nuclear matrix element 
Mk:- 

To first order in the Fermi momenta, the energy difference 



AE = Eo- ^int = AEo + k-W + q-V 



Here 



AEo = Erjipr,) + Ea{-Pri) - Eh{-\pn) - En{\pd) - E^{\p^ - \pd) 
and the relativistic relative velocity vectors V and W , 
V = 



(15) 



(16) 



\pd) 



4K(|a, 



■Pv 



WPd) 

1 
2 



+ 



1 



W 



EAjPri - \Pd) Er,{\pd) 
-^A\Pv - \Pd) + ^h{-\p^) 



Pd 



EniiPv - \Pd) Ehi 



iPv) 



Pr, + 



2E^{%-^^Pd) 



Pd , (17) 



depend only upon external kinematic variables. 

After replacing rrin/ En{\pd) by rrid/ Ed{pd), the linearisation simplifies the 
matrix element A4 of eq.(|l) to 



M 



frih 



{2TxY2E^{lPr,-\pd) Ed{\Pd) Eh{-\pr,^ 



NT\f^m~m'd 



X 



j d^k J d\ 



AEo + k-W + q-V 



M 



K 



le 



;i8) 



Close to threshold the internal helion is almost at rest but relativistic effects 
are important for the incoming deuteron. Taking them into account leads to 



M 



1 



1 



(2T:f2E^{{p, 



-jPd, 



Nt \f^rn~m'd j d^k j d^Q 



-M 



AEo + k-W + q- V + ie 

(19) 



K, 
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where 



v; 



(20) 



For production in the forward direction, the vectors V and W are parallel to 
the beam direction in the overall cm. frame, 



V 



-pd, W = pd. 



(21) 



After neglecting the Fermi momentum, the pd 7r° ^He amplitude should be 
180°. 

To evaluate the matrix element Ai of eq. (|T^) , we first rewrite the denominator 
as an integral over a parameter t, 



taken at 



M 



(27r)3 2E^(|p,-lp,) 



Nt V4:m^md dt J d^k J d^q e''^^^+''^ Mk 



(22) 

which allows the integral over the Fermi momenta in the wave functions to be 
evaluated analytically. 

The spin-algebra can be handled by the same tensor techniques used in I. For 
collinear kinematics, (0), the matrix element yields an expression for the F 
amplitude of eq.(0) in terms of a single complex form factor 



5(H>,y') = (27r)ViV^y^ dte''^''y:i-tW)l^^s{tV') + -j=ipDitV')j . (23) 

This involves an integral over configuration-space ^He and S- and D-state deuteron 
wave functions. 

The F amplitude may be written in terms of this form factor as 



NtGBS 



(24) 



This is the expected structure since both the B and F amplitudes correspond to 
incident deuterons with helicity m = ±1. 

At threshold, we make use of eqs.(|l^) and (|12|) to replace the squares of 
the amplitudes Q and B by the corresponding (polarised) cross sections. From 
eqs. (pi) (0) it follows that 



Pd da 4 

— -j7z{dd^ Her)) 

Pr, d\l 



Pd da 
Pt, d\l 



64 



3 {27.YEiapd) 



YM\S{Wy)\^x 



-1 m=±l r 



p^ da 

— :T7T(7r p^r]n) 

Prj d\l 



(25) 
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where the value Nt = 12 has been substituted. 
The mass factor Ym is defined as 



— 5 7 \ l^Oj 

The two factors of the expression for Ym are related to the cm. energies of 
the participating vr+n rjp and dp — » ^He 7r° reactions, and originate from terms 
in the cross section formulas ([l^ and ([I^). It will be noted that the approximate 



threshold formula is similar in structure to the corresponding one for the 
near-threshold pd ^Herj reaction derived in I. 

In eq. (|25| ) the value of the dp ^He 7r° cross section should be taken at 
Opj^ = 180", for deuteron helicity m = ±1. Such cross sections have been well 
measured in our energy domain fl^ . 
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4 Final state interaction 



Though the experimental amphtude for pd ^Herj varies very rapidly near 
threshold |Q, the predictions of our two-step model are much smoother since the 
form factors hardly change for energy steps of a few MeV. This sharp dependence 
was ascribed to a strong final state interaction and in I we developed a 
model for the S-wave enhancement factor Q based on an effective wave technique. 
Generalising this to the case of four nucleons, it is straightforward to see that the 
effective wave "^ii^Xi) incident on a nucleon at Xi is given by 



for which the denominator is 



AM 



1 + / — ^gifc-(x,-x,) _ 2 



D 



(27) 



D= 1 + f 



I 



1-3/ 



ikl 



I 



(28) 



In the above, / is the rjN scattering amplitude and k the momentum of the rj 
scattering off the nucleons in the ^ife-nucleus, which is assumed to be represented 
by four scattering centres placed at the vertices of a tetrahedron with side length 
/. 

The point '^He charge form factor in the rigid tetrahedron model is 



^ch(g) = Jo qi 




(29) 



To obtain consistency between the '^He r.m.s. radius and the form factor minimum 
[ P^m , the Bessel function may be smeared over the interparticle separation with a 
weight function 



G{1) = N{1- Vn) 9(1 - l^^) expi-\\l - Af) , 
with parameter values 

Wn = 0-50 fm , A = 0.82 f m , A = 0.543 fm"^ . 



(30) 



(31) 



The subsidiary maximum in the form factor is then underpredicted by about a 
factor of 1.7, but this feature is sensitive to meson exchange- current effects. 
The jS-wave enhancement factor becomes 



D 



Akl 



1 + / — f,ikixi-xi) _ 2 



(32) 
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Averaging over the orientation of the tetrahedron projects out the S-wave, to 
leave an expectation value with respect to the lengths of its sides 



= 



dlG{l) -h+/_{3jo(M)-2} 



) 



(33) 



Following the prescription used in I, for finite nucleon masses we take the 
effective amplitude 



frfN is here the 77-nucleon cm. scattering amplitude, evaluated at a cm. momen- 
tum, 



where k is the cm. momentum in the '^Her] system 

In fig. 2 we show the variation of | f2 p withp^ in the illustrative case where the 
77-nucleon scattering amplitude is taken to be constant at the scattering length 
value a — (0.476 -I- 0.279i) fm used in I. At threshold the attractive r^-nucleon 
interaction enhances the predictions by a factor | VL p= 1.9. However the rapid 
fall-off from enhancement to suppression around 25 MeV/c is very striking. 




(34) 




(35) 
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5 The n ^He ^He 77 reaction 

If the near-threshold energy dependence is characteristic of an ri'^He final state in- 
teraction, then the effect should be largely independent of the particular entrance 
channel. An alternative means of reaching this final state is through studying 
the n ^He — >■ '^He i] reaction or the equivalent one with a proton beam and tritium 
target. 

Just as for the dd case, there is only one amplitude at threshold or in the 
forward direction, and it may be written as 



In terms of the amplitude Ai, the differential cross section becomes 



^{n^He^ 'Her,) 



The two-step model for the n^He '^Herj reaction shown in fig. 3, merely 
involves a rearrangement of the initial nucleons in fig. 1. Thus M.k contains the 
same type of matrix elements as eq.(|^) with the exception of the initial wave 
function, 

-Mk = E ^i^ic ^ p'J)Mi7i% ^ p,rif)M{p,t ^ af)Minid tiH) . (38) 

int 

After carrying out the linearising procedure as in §3, it is straightforward to 
show that the two-step contribution to the amplitude H is 

H = J^.t^, {SsiW, V') - V2SniW, V'))+^^ {Ss{W, V') + Sn{W, V')/^) 

(39) 

where an extra factor of three arises from including also the analogous charged 
pion contribution in fig. 3. 

The form factors are defined by 

Ss,d{W, V ) = (27r)3 ^N^N, e^*^^«C(-^W^)05,D(t^ ) , (40) 

where 05 (r) and 4'd{''^ are the helion S- and D-state wave functions of the pd 
partition, a parametrisation of which is given in ref . fl^ , and the number of such 
partitions A^^^ ~ 1.5. 

Here the bare intermediate state energy difference is given by 

AEq = E^(p^) + Eo,{-p^) - Et{-\p^) - Ep{-\pn) - E^Hp^ + \pn) (41) 
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and the relativistic velocity vectors W and V are defined as 



V 
W 



(42) 



A crucial difference with the dd '^Herj calculation is that the input pd 
■^Hen^ amplitudes are required in the forward direction but at a rather higher 
incident energy. The experimental data |]l5i indicate that 



A{9p^ = 0) p » I B{9p^ = 0) 



(43) 



If we only keep the dominant A amplitude, then we can derive a linear relation 
between the n ^He — *• ^He rj and the unpolarised input cross sections 



Pn do_^ 

Prf dfl 



n He 



128 {2n)^Elilpr. 
Pd da 



Ptt dQ 



{pd 



Yi,\Ss{W,V')-V2SD{W,V')\' 
Pn da 



Prj dVt 



{tx p rjn) 



The mass factor Y^.^ is similar to that given in eq. 



Y' 



{nip + m-riY ml + pm^(2m„ + rUr,) 



{ma + m^)2 



(45) 



It is of course to be expected that the final state enhancement factor | f2 P 
should be the same here as for the dd ^He rj reaction. 
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6 Predictions and Conclusions 



The Breit-Wigner fit to the it p ^ rjn threshold amphtude given in 
Ptt da 



dfl 



{it n rjp) 



(680±60)^6/sr 



yields 
(46) 



where the error bar is taken from the direct measurement of Binnie et al. . 

At the dd ^He t] reaction threshold, the laboratory deuteron kinetic energy 
is 1.12 GeV, and for this energy the helicity m = 1 differential cross section of 



the pd 



He 71 reaction at 9] 



pd 



180 is measured to be lib 



Pd da 
Pt, dfl 



{pd ^He 7r°) 



M=l 



;i.0±0.1)/ib/sr 



(47) 



Using the Gaussian parametrisation (|^) of the ^He : n^He cluster wave func- 



tion, together with the Paris deuteron wave function [|1^, the magnitude of the 
transition form factor of eq. (P3| ) is found to be 

(48) 



\S\ = 11.2 fm"' 

at the threshold for //-production. 

Before introducing final state interactions into our model, eq. 
threshold amplitude squared of 



predicts a 



f{dd^ ^HeT]) 



Pd da 
— -77z{dd 
Prj a\l 



(44±8)nb/sr, 



(49) 



to be compared with the experimental value of about 30 nb/sr. 

The energy dependence of the two-step form factor and the input dp —>■ ^He 7r° 
cross section are both negligible for small Prj. The only significant variations arise 
from the final state interactions and, to a lesser extent, the ii'p rjn amplitude 
squared. 

After including the final state interaction effects, using the formalism of §4, the 
predictions shown in fig. 4, which for ease of comparison have been multiplied by a 
scale factor of 0.75, seem to have too sharp a momentum dependence as compared 
to the SPESIV data 0. Though these dat a were taken over a very small range 
in prj, and the lowest point is subject to large corrections due to energy losses in 
the target and similar effects, the much smaller slope as compared to pd —>■ ^He rj 
is confirmed in the SPESIII experiment over a much wider domain in p.^ 0. It 
should however be remarked that the FSI factor depends sensitively upon 
the //-nucleon scattering length and it has recently been suggested that it 
might be even larger than that which we have taken here. 

In the case of the n^He — > '^Herj reaction, the neutron laboratory kinetic 
energy at threshold is 0.75 GeV so that the dp — > ^ii/evr^ input is required at 
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laboratory kinetic energy of = 1.50 GeV for 6*^^ = 0°. Using the data of 
ref.[jl5|, we find 



Pdda 



- unpol 



(2.8 ± 0.3) ;ub/sr . 



(50) 



Taking the Germond- Wilkin parametrisation of the ^He : dp cluster wave 
function |T^, we find the magnitude of the form factor to be 



Ss-V2S_ 



D 



17.3 fm^ 



so that the predicted amplitude squared becomes 



Pd do_^ 
Pri dfl 



n^He 



(51) 



(21 ± 4) nb/sr , (52) 



which is about a factor of two smaller than the dd prediction. One of the major 
reasons for the difference between the two predictions is the much larger numerical 
factor in eq.(p5D compared to (^). In the dd case there are four times as many 
graphs, all of which contribute coherently. 

Unlike the dd case, the two-step model for the n^He '^Het] reaction in fig. 3 
should have contributions from intermediate spin-singlet NN (d*) states. From 
the experience gathered in I, these could well increase the prediction by over a 
factor of two, to make it comparable to the dd result. Nevertheless, it is much 
harder to investigate this reaction experimentally due to the difficulties of working 
with a broad-band neutron beam or a radioactive tritium target. Attempts at 
such a measurement at LAMPF were aborted due to problems with the tritium 
target 

There are however few other possibilities. The low energy r] ^He system could 
be investigated through pion production or even more cleanly through the photo- 
production 7 ^He —>■ 1] ^He since sufficient energy resolution can now be obtained 
P^ . On the other hand photoproduction of the rj'^He system is strongly sup- 



pressed due to the dominance of isovector photon coupling to the A^*(1535). 

It should be noted that the near-threshold kinematics are not as 'magic' 
as they are for pd ^Hei], since in the case of dd ^Her] the zero- Fermi- 
momentum energy difference Ai^o = —135 MeV while for the n^He — > '^Herj 
case, AE'o = -1-140 MeV. The imaginary parts of the form factors become rather 
important, and the situation is quite similar to the pd —>■ ^He [uj, rj') cases, which 
are successfully reproduced by the two-step model 

In contrast to the pd — ^ ^Herj estimate, our calculation of the dd ^Herj 
amplitude is slightly higher than experiment. Since contributions from all con- 
stituents in the nuclear wave function have been included, this reinforces the 
suggestion made in I that the pd discrepancy is due to the omission of the inter- 
mediate d* effects. Of course initial-state damping, which we have not considered, 
could well reduce somewhat the ^ife rj production rate. 
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Figure Captions 



Figure 1 : Double-scattering diagram for dd ^He 1] in terms of the constituent 
pd ^Hen^ and ir^n rjn amplitudes. Single lines are nucleons, double 
deuterons, triple ^He and quadruple "^He, the dashed meson lines being labelled 
individually. The momentum labels refer to the overall cm. system. In total 
there are eight such diagrams corresponding to the interchanges n, ^He — ^ p, ^H, 
7r° — > 7r+ and exchanging beam and target deuterons. 

Figure 2 : Enhancement factor | f2 p of eq. (|33[) evaluated in the constant ?7-N 
scattering length approximation. 

Figure 3 : Double-scattering diagram for n'^He —>■ "^Herj in terms of the con- 
stituent pd ^He and ir^n — > rjn amplitudes. 



Figure 4 : Spin- averaged amplitude squared |/|" 



Pri d\l 



for 



the dd — > ^Herj reaction near threshold, as a function of the rj cm. momentum 
Pri- The theoretical curve of eq.(P), with FSI enhancement and multiplied by a 
constant factor of 0.75, is compared with experimental data taken from ref. H]. 
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Fig. 3 
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